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Abstract - 

A self consistent field theory for compressible polymer mixtures is devel- 
oped by introducing elements of classical density functional theory into the 
framework of the Helfand theory. It is then applied to study free surfaces of 
binary (A,B) polymer blends. Density profiles in the one- and two-phase re- 
gion are calculated as well as chain end distributions and chain orientations of 
the minority and the majority component. In the ideally symmetric mixture, 
in which the individual properties of polymers A and B are the same and both 
have the same surface energy, polymers of the minority component segregate 
to the surface, where they are exposed to less polymers of the majority com- 
ponent. This effect can only be captured correctly, if one accounts for the fact 
that the monomer- monomer interaction has finite range. As a consequence, 
the Flory-Huggins-parameter varies in space and depends on the concentration 
profiles and their derivatives. The surface segregation calculated with such an 
ansatz, without any fit parameter, is in reasonable quantitative agreement with 
data from recent Monte Carlo simulations. 
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1 Introduction 



Many important properties of materials are determined by the structure of 
their surfaces (adhesive properties, thin film properties etc.). Polymeric sub- 
stances are affected by the presence of a surface on many different length scales: 
rearrangement of monomers similar to the layering and packing of molecules 
at surfaces of simple liquids, reorientation of monomers, reorientation of whole 
chains, surface induced ordering in copolymer systems, wetting and other sur- 
face induced transitions [[l], 0. One of the most remarkable effects at surfaces of 
polymeric mixtures is the phenomenon of surface segregation - the composition 
of the melt close to the surface differs from that in the bulk. The reasons for 
such a behavior are manifold. Shorter chains tend to segregate to the surface 
because they loose the least entropy entropy there 0, |, ||; Stiffer chains are 
favored at the surface because they pack more easily ||. Hard walls may inter- 
act with the different polymer species with different interaction strength. At 
free surfaces, segregation is induced by "missing neighbor effects" : The poly- 
mer species with the lower monomer-monomer attraction segregates, since it 
is less penalized for the lack of neighbors at the surface - it's surface energy 
is lower 0. A similar effect occurs even in completely symmetric mixtures of 
two incompatible, but otherwise "identical" polymers. The minority compo- 
nent enriches at the surface, because it has less contacts with polymers of the 
majority component there. This interaction driven segregation is a general 
phenomenon, which is found in many systems with a tendency of bulk order- 
ing (the "ordering" here being the demixing of the two polymer species). It 
can be expected in the one-phase region as well as in the two-phase region, 
and is strongest in the vicinity of the critical demixing point 0. An additional 
aspect of surface segregation in polymeric systems has to do with the competi- 
tion of length scales there: the single chain gyration radius competes with the 
typical length scale for collective concentration fluctuations, which diverges at 
the critical point. 

Experimental studies have mainly been concerned with surface segregation 



in polymer mixtures in the one-phase region J^-[T1J. Popular systems are, 
for example, mixtures of polystyrene (PS) and deuterated polystyrene (dPS), 
or of hydrogenated and deuterated polymethylmethacrylate (PMMA and dP- 
MMA). The surface of PMMA/dPMMA mixtures saturates with dPMMA 
almost independently of chain lengths and bulk composition JO]]. Hence the 
difference between the surface energies of the components is presumably very 
high and entirely dominates the segregation process. The segregation tendency 
of dPS in PS/dPS mixtures is much weaker and the dPS excess at the sur- 
face depends strongly on the bulk volume fraction. The relative importance of 
different segregation mechanisms can be estimated by simple considerations. 
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The excess surface energy of dPS at vanishing bulk volume fraction is roughly 
(i pd 2 — 3 ■ lCr 3 nm|| at temperature 184°C The Flory Huggins parameter 
X, which measures the relative repulsion of monomers of the two components, 
takes the value x — 1-5 • 10~ 4 , and the polymer segment length is a = O.Q7nm. 
From comparing \x to x a i one concludes that the segregation process is still 
mainly driven by the excess surface energy, but that interaction driven seg- 
regation may become important. An even stronger effect can be expected in 
mixtures where x is larger. 

Apart from experiments, Computer simulations provide additional under- 
standing and a good testing ground for theories |12fl. Monte Carlo simulations 
of simple one component polymer melts at surfaces have given insight into the 
structure of polymer chains close to the surface - chain end distributions, chain 
orientations, etc. ||13|| -||16||. Being computationally very demanding, simulations 
of polymer blends at surfaces are still scarce. Wang et al have investigated the 
wetting behavior of very short chains |17]] at walls with very strong attraction 
of one species. In a recent study, Rouault et al[18j consider incompatible sym- 
metric polymer blends at "neutral" walls and find, as expected, segregation of 
the minority component there. In contrast, Cifra et al|19[] study fully compat- 
ible blends, where the missing neighbor effect leads to the enrichment of the 
majority component at the surface. 

Theoretical treatments of polymeric surfaces have often dealt with simple 
Flory-Huggins-de Gennes functionals of the concentration profiles in "square 
gradient" approximation]^]- [Q. Jones et alJITJ, first pointed out that 
this type of theory fails to reproduce the correct form of segregation profiles: 
Unlike profiles typically obtained from Ginzburg-Landau theories, experimen- 
tal profiles are flattened right at the surface. More sophisticated mean field 
approaches take the chain connectivity into account more explicitly, either in 
the framework of integral equation theories or density functional theories [f25|1 - 



33fl , or in the self consistent field theory originally developed by Helfand and 
othersf34|-|3^], or, lately, in combinations of the two approaches [|3l], |33f . A 
recent study by Genzer et al|39[] has shown that profiles of polymer blend sur- 
faces calculated within the self consistent field theory can indeed be fitted to 
experimentally measured profiles in a satisfactory way. However, Genzer et al 
treat the blend as an incompressible mixture, a simplification which is clearly 
not justified when treating free surfaces. Surface segregation is enforced by 
external ad hoc potentials acting on the surface layer. The theory does not 
allow for a microscopic treatment of missing neighbor effects. 

The importance of compressibility effects in polymer melts has been recog- 
nized by several authors. In the early work of Helfand et al, they are taken 
into account within a quadratic approximation. This is reasonable as long 
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as one deals with small total density fluctuations, not in the vicinity of a 
surface where the density goes down to zero. Some work has been devoted 
to elaborated density functional theories of one component polymer melts at 
surfaces p5|, |26| , [27| , |32| . Those studies have had remarkable success in repro- 
ducing the density profiles measured in corresponding Monte Carlo Simula- 
tions; however, they mainly focussed on the microscopic structure of a melt 
over the length scale of the monomer size - here, we are interested in composi- 
tion variations on a larger scale. Freed and coworkers have developed density 
functional theories of compressible polymer blends |3(J. Recently, Hariharan 
and Harris have presented a self consistent field study of the surface of a com- 
pressible copolymer melt within the framework of the Scheutjens and Fleer 
lattice theory |37|, |40fl . Compressibility can be introduced into lattice models in 
a relatively straightforward way by admitting vacancies as a new type of "par- 
ticles" which may occupy a lattice site. To the present author's knowledge, 
no attempt has been made so far to study surfaces of compressible polymer 
mixtures in continuum space. 

This is the scope of the present work. On the base of the equation of state 
which is assumed to be known and using elements of simple density functional 
theory, a self consistent field theory for compressible polymer melts will be 
developed. In self consistent field theory, monomer-monomer interactions are 
commonly treated as contact interactions V(r) oc x^if)- Such a potential 
Ansatz is not suitable for the description of important missing neighbor effects 
- e.g. the effect that, close to a hard wall, monomers a z > have no neighbors 
at z < 0. Hence the theory will be extended such that it accounts for the finite 
range of interactions. The parameters of the theory will be adapted to the 
bond-fluctuation model used in the Monte Carlo simulations of Rouault et 
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The bulk properties of the model are very well known from previous work[[4l"|, 
42f . Miiller et al performed large scale simulation studies of interfacial prop- 



erties in incompatible polymer mixtures^]. The comparison of his results to 
self consistent field calculations showed that the theory is rather successful in 
predicting many quantities even at the relatively short chain length of N = 32, 



except for the interfacial width[]45|. The discrepancies between theory and sim- 
ulation becomes smaller as the chain length is increased, but the convergence 
is extremely slow. Similarly, we will see here that our self consistent field 
calculation yield the quantitatively correct values of the surface segregation, 
without any fit parameter, but fail to predict the correct decay lengths of the 
profiles. 

The paper is organized as follows. The theory is developed and the model 
parameters are determined in the next section. In section 3, the surface seg- 
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regation at neutral walls is calculated in the one-phase and in the two-phase 
region and compared to the available simulation data. The influence of chain 
length asymmetry on the surface segregation is also studied. Then the struc- 
ture of the blend at the surface is analyzed in more detail. The evolution of 
density profiles as one moves away from the critical point will be examined, 
chain end distributions and chain orientations will be calculated. The results 
are summarized in section 4. 



2 Self consistent field theory for compressible melts 



We consider mixtures of flexible polymers A and B in the vicinity of a surface. 
Each chain of type i (i = A or B) contains monomers and is characterized 
by a statistical segment length 6 i; The radius of gyration of a chain in the ho- 
mogeneous melt is thus given by R g ^ 2 = 6j 2 A^/6 (random walk). Unless stated 
otherwise, we will assume that polymers A and B have the same properties, 
i.e. N, b and the bulk density in a pure system p b are equal for both species. 

Molecules are treated as space curves r(s), with s varying from to 1. The 
partition function for a system of ha polymers of type A and % polymers of 
type B has the general formp5 



z =^-, / n^(-)}n^(-)}exp[-^{p A ,M], 



1=1 j=l 

where j3 = is the Boltzmann factor and pa, Pb are monomer density 

operators 

Pa.b(?) = N Xf t ds 5(f- n(s)) (2) 

and ^{pajPb} represents a coarse grained free energy functional, which will 
be discussed in more detail below. Individual space curves in the functional 
integral are assigned a statistical weight V{f(-)} = V{f(-)} Vg{t{-)} with 

^{r(-)}=A/-exp[--4- / ds\U% (3) 



2Nb 2 Jo 'ds 1 

with the normalization factor N '. Flexible chains are hence modelled as con- 
tinuous space curves, which are distributed according to a Gaussian stretching 
energy with spring constant 3/Nb 2 . 

We treat the problem in mean field approximation by replacing the monomer 
density operator p(f) in eqn (0) with the average monomer density Pi(r) := (pi(r)). 
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Monomer-monomer correlations are ignored. The chains can then be treated 



as independent random walks in the external field Wi{r) = [ |35| |. It is 



useful to define the end-segment distribution functions 

Qi(f , s)= f V{f(-)} exp[- f ds'WiirWMra - r(s)) (4) 



which obey the diffusion equation [p4| 

{ k^s-\ h2 ^ +Wi}Q ^ s) = {} (5) 

with initial condition Qi(r, 0) = 1. The function Qi gives the statistical weight 
for a part of a polymer of type i and length sN, with one end fixed at position 
r. The average density of type i monomers (with orientation u) is then given 
by 

Pi(r) = Pb ds Qi(f, s)Qi(r, 1-s), (6) 



which completes the cycle of self-consistent equations. The calculated profiles 
can be used to calculate the configurational entropy of a polymer of type % 

S t = knilogZi + - [ drWi{r) Pi {r)), (7) 

Hi J 

where the single chain partition function Zi — J dfQi(r,l) = Nrii/pb is a 
constant. From there one gets the free energy 

@F = ^T{p a ,pb) - -r- {n A SA + n B S B ). (8) 

k B 

We now turn to the discussion of the free energy functional JFjp^, p^}, 
which defines the concrete model. We assume that species dependent interac- 
tions between monomers can be treated as perturbations of a reference system 
"0" , in which there is no distinction between monomers of different type. 

T{p A , Pb] = Fo{p} + Ti n ter{pA, Pb} with p = Pa + Pb (9) 

More specifically, our Ansatz is based on a reference system of noninteracting 
polymers with simply hard core repulsion. We can treat the reference system 
in local density approximation 

{3F {p} = ( drp{r){f[p{r)} - p). (10) 



The bulk free energy per monomer f(p) can be calculated from the equation 
of state n(p) (with the pressure II) using 

f(p)= f P dxU(x)/x 2 - (11) 
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the chemical potential p drives the total number of particles and is related to 
the bulk density p& via p = f(pb) + Pb§^- We will use the equation of state 
resulting from the generalized Flory theory in continuous space proposed by 
Dickman and coworkers f4p|, 4"7[ . 



n = - 1], (12) 

for long enough chain lengths. The constant Cq is Cq = 0.6583, and the 
"packing fraction" rj is given by rj = a 3 p7r/6 with the "monomer diameter" 
a = 1.96oq in the bond fluctuation model [^7] (where a is the lattice spacing). 
This theory has been remarkably successful in predicting the correct pressures 
for the athermal bond-fluctuation model over a wide range of densities. At the 
interesting density p = l/16ag 3 , however, it yields a compressibility k = 34a 3 , 
which is slightly too high (k « 31a 3 ) according to simulations! 42]). In order to 
correct for this, a higher order term is added to f(p), 

/M = C " (4 " 3 "^" 3) with 6 = 3.0. (13) 

(i -vr 

The ad-hoc correction to the generalized Flory theory (b = 0) is motivated 
by the fact that the equation of state derived from (|T3"D, II = p 2 f'(p), is 
numerically almost identical with ( |i~2"|) in the relevant density region < p < 
Oq/16. Calculations have also been performed at b = or with a completely 
different form of f(p) (a quadratic expansion around p = Oq/16, cf. |4"5| ) - the 



results appear to be rather insensitive to details of the equation of state. 

The formalism so far takes care only of the hard core repulsion between 
monomers. In a perturbative treatment, monomers A and B are assumed to 
interact with each other with additional integrable potentials Uaa{^)i Ubb{R) 
and Uab{^), and the energetic contribution to the free energy is given by 



(14) 



where the indices i,j run over the polymer species A and B, and p[f(r, r") is 
the pair density of type % monomers at point r and type j monomers at point 
f . We use the simple approximation 

p[f(fy) = j(f- r*)pi(?)Pj(f), (15) 

where 7(f) does not depend on densities any more, but does account for local 
effects of chain connectivity. A possible choice is 7(f) = inside the excluded 
volume of the central monomer, and j(f) = 1 — p(f) outside, where p(f) is 
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the probability that direct neighbor monomers along the chain block point r 
for occupation. The Ansatz ( ITS] ) is chosen in the usual spirit of Flory Hug- 
gins type theories, where monomers are subject to averaged fields created by 
all other monomers, and the effect of long range concentration correlations 
induced by chain connectivity is ignored. This implies that the demixing be- 
havior is governed by the interchain interactions, as has indeed been found in 
simulations, see ref. fP]] . Hence the pair correlation function 7(r) is taken to 



be independent of the identity of the monomers. Note that, as soon as intra- 
chain interactions become important, demixing goes along with the collapse 
of polymers and the Helfand treatment of chains as random walks is no longer 
appropriate. 

Equation ( |14|) can now be rewritten as 



Pinter = J dfdr'^f - T*)\p{f) p{f)V x {f - r 1 ) + p A (r)p B (f^) V 2 (f - f) 

+ {pa{t) - PB(r))p(r)V 3 (r - r)} (16) 

with 

Vi(r) = ^(U AA (r) + U BB (r)) 

V 2 (r) = ^(2U AB (f) - U AA (r) - U BB (r)) 

Vs(f) = ^(U AA (r)-U BB (r)). 

The first term depends only on the total density profile p(f) and does not 
introduce qualitative changes compared to the reference state JF . The second 
term describes effective interactions between monomers A and B, and the third 
term an effective field favoring one of the components. One can see this from 
looking at the homogeneous bulk system. For i = 1, 2, 3 we define 



Ub,i = Pb J df-y{f)Vi{f). (17) 
The contribution (|16"D to the bulk free energy is then given by 

(3Finter,b = «b,l n + M&,2™4™b/™ + U h ^{n A - U B ), (18) 

with the total number of monomers n A , n B and n = n A + n B . Hence we can 
identify Ub, 2 with the Flory Huggins parameter Xb, and Ub,3 with an effective 
chemical potential difference hb- 

We assume that the potentials Uij are short range and expand the profiles 
PA,B(f) around r. Since we study a planar surface, the system is inhomo- 
geneous only in one space variable z. Using / dfj(r)Vi(r)z = (due to the 
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inversion symmetry of V and 7) , and defining 



k i = 7 p- I df 7 (f)VHr> 2 , (19) 

one can reexpress terms of the form J dfdf^ p a (f)pj3(^)Vi(f — r > ) , y(f ~ f') by 
J dzp a (z)(pp(z) + p'p(z)ki)ub t i/ pb, where = A or 5. One obtains the total 
free energy p9| 



+x(z)PA{z)p B (z)/p b + h(z)(p A (z) - ps(^))} (20) 

with the space dependent Flory-Huggins parameter x( z ) = X&(l+2^2(^+^)) 
and the space dependent field 7t(z) = /i&(p + k 3 p")/pb. The difference to the 
bulk values of x and h reflects the missing neighbor effects and is ultimately 
responsible for the occurrence of segregation. This is obvious for the field h(z) 
- as we will see, it is also true for x( z )- From eqn (PD| ) one can calculate the 
mean fields 

W A {z) = ((z) + Xb(pB + k 2 p'^/ Pb + h(z) (21) 
W B (z) = ((z) + X b(pA + hp A )/pb-h{z) 

with 

df 

C(z) ~pf p + f(P) -» + Cinter{z), (22) 

(inter(z) = 2(p + fcip")— + ((pA ~ Pb) + h(p A ~ Pb)) — 

Pb Pb 

We will assume that the perturbation Jointer is small compared to JF and that 
the corresponding contributions to ((z), Qnter(z), can be neglected. 

To summarize, elements of a simple classical density functional theory (the 
local density approximation for the treatment of hard core monomer inter- 
actions in connection with perturbative treatment of attractive interactions) 
have been introduced into the Helfand formalism, in order to enable the study 
of compressible systems with intermonomeric interactions of finite range. A 
similar approach has been taken by S.K. Nath et al in order to study surfaces 



of homopolymer meltsp3|. We note that it has some inconsistencies, polymers 
are treated as gaussian strings, whereas monomers are assumed to have finite 
size. This should however be of little importance, as long as we are interested 
in density variations on scales larger than the monomer size, and as long as 
the chains are much longer than the monomers. 
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The model parameters are adjusted to fit the bond fluctuation model, with 
the model parameters used by Rouault et al in their Monte Carlo simulation |18[|. 
The bond fluctuation model is a lattice model, where monomers occupy each 
a cube of eight neighboring sites and are connected by bond vectors of length 
2 < d < ylO lattice spacings. Thus the excluded volume of a monomer 
is V exc i = 27 Oq. Simulations are often carried out at monomer density p = 
l/16ag 3 , where one half of all available lattice sites are occupied and the model 
reproduces many properties of dense polymer melts. Rouault et al simulated 
polymers of length 32 and used the interaction potentials Uaa{t) = Ubb(j) = 
—Uab{^) = —ksTe for \r\ < a/6, i.e. they have V\ = — e/2, V 2 = 2e and 
V3 = within the interaction range. No effective field h(z) is induced, since 
the interaction parameters are symmetric with respect to A and B. The sites 
in the interaction region 2 < \r\ < y/6 of a monomer are at the relative posi- 
tions r — (200), (210), (211), plus permutations of the coordinate axes and sign 
reversals. Hence the interaction region covers N t = 54 sites, A2 = 18 of which 
are located at distances z = ±2 from the central site along the z-direction, 
and Ni = 24 at distances z = ±1. Assuming that a site in the interaction 
range of a monomer is blocked by one of its neighbor monomers along the 
chain with the probability p = 1/3, a rough estimate of the function 7(f) 
yields 7(f) = 1 — p = 2/3 in the interaction region and therefore Xb = ^2^// 
with the effective coordination number z e ff = (1 — p)N t pb with p& = 1/16 (cf. 
eqn.flTTI)). The value z e ff = 2.25 calculated from this approximation is reason- 
ably close to the average number of interchain contacts found in simulations, 
z e ff = 2.65 at chain length 32[^]. Similarly, one estimates from eqn ( |T?D and 
eqnflg) k 2 = £j fN j /(2N t ) = 8/9. 

It should be emphasized that the theory was adjusted to the bond fluctua- 
tion model, a lattice model, for the sole reason that simulation data are avail- 
able, which permit to test the theoretical predictions. Adjusting the theory to 
continuous off-lattice models causes no conceptual difficulties in most cases. It 
requires the knowledge of the equation of state in a reference system of identi- 
cal polymers, a reasonable guess of the pair correlation function j(r), and the 
knowledge or a good ansatz for additional integrable intermonomeric interac- 
tions Uaai Uab and Ubb, which are small enough to be treated perturbatively. 
Extensions of the theory will be necessary if the origin of the Flory-Huggins 
parameter \ is mostly entropic, e.g. if polymers A and B demix as a result of 
different monomer sizes. 

A surface at z = is imposed by requiring Pi(z) = for z < 0; at the bulk 
side, the boundary conditions ^Qi{z max ,s) = and Wi{z max ) = W i)bu i k were 
used, where Wi^uik denotes the bulk value of the fields Wj. The numerical 
method used to find the self consistent solution of the diffusion equation was 
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an iterative relaxation technique described in ref. P5| . The iteration cycle 
was stopped after reaching an accuracy r = J2% I dy(AWi) 2 of 10~ 14 . (relative 
accuracy > 1CT 15 ). 



3 Results 

The system shows the usual bulk demixing phase behavior. From eqns (|5]) and 
d||), one can see that the bulk volume fractions $a,b = Pa,b/ Pb obey the mean 
field equations $j = exp(— NiWi) or, equivalently, dF/dpi = log($j) / 'N{+Wi = 
with the Flory type bulk free energy per volume 

(3F = ( PA + PB )f{ PA + PB ) + Y,%r M^) + —papb - E H 0i> (23) 

i Pb i 

and H A = p + h b + l/N A , H B = p — h b + 1/N B . It is useful to rewrite the 
free energy as a function of the total density p = p A + Pb an d the density 
difference m = p A — p B - A homogeneous blend is stable or metastable with 
respect to small fluctuations in m as long as [d 2 F/ \dm) 2 } p = (Xs~x)/(^Pb) > 
with the spinodal \s = 0-/{Na$a) + 1/ (N B § B ))/2. The spinodal ends in the 
critical point, which is characterized by [d 2 F/(dm) 2 ] p = and [d 3 F/(dm) 3 ] p = 
((1/(N a $a) ~ l/(^$|))/(2pb) = 0. Hence one gets the critical \ parameter 
and critical volume fractions 




(24) 

At low temperatures or x > Xc, the mean field equations have no stable ho- 
mogeneous solution for a range of concentration differences m, i.e. the system 
demixes. This implies that there exist two distinct solutions in a range of hb, 
one finds metastable states or two phase coexistence; at Xc the two solutions 
merge continuously into one. The order parameter of the system can be de- 
fined as the difference of volume fractions in the two coexisting phases. In the 
symmetric system with equal chain lengths A^^ = iV^, two-phase coexistence 
occurs below the critical point (x-/V) c = 2 at hb — (see Fig. 1). 

3.1 Surface segregation 

Figure 2 shows a segregation profile in a symmetric mixture at bulk two phase 
coexistence. The volume fraction of the minority component is increased al- 
most by a factor of two at the surface. After a flatter part, the composition 
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profile decays towards the bulk value within 2 — 3 times the gyration radius 
R g . The total density profile is shown in the inset: The density reaches the 
bulk value much more quickly than the composition profile. This reflects the 
different length scales involved in the problem: The relaxation of the total 
density profile involves rearrangement of single " monomers" , it is hence driven 
by the Kuhn length b. The composition change at the surface, on the other 
hand, goes along with the segregation of whole chains - the thickness of the 
segregation layer has to be of the order of at least the gyration radius R g . The 
fine structure of the profiles will be discussed in more detail later. It is worth 
noting that the extent of surface segregation strongly depends on the range of 
the monomer interactions or the value of In the limit of k 2 = or pure 
contact interactions, segregation is almost entirely suppressed (dashed line). 

Apart from the bulk phase diagram, Figure 1 also indicates the surface vol- 
ume fraction of component A in the different phases as a function of "temper- 
ature" 1/xN a t two-phase coexistence. The difference between surface 
and bulk volume fraction predicted by the theory is compared to the simulation 
results of Rouault et al in Fig. 3. (Note that our mean field theory overesti- 
mates the critical temperature and does not describe correctly the Ising-type 



critical behavior found in short chain systems PT|, but the region near crit- 
icality is not our main interest here.) The simulations were performed in a 
slab geometry with a slab thickness of around three times the gyration radius. 
In such a thin film, the critical temperature is reduced, and so is the order 
parameter in the two phase region. The difference between the surface volume 
fraction of the minority phase in the thin film and the bulk volume fraction in 
an "infinite" system can hence be taken as an upper bound for A$^. A lower 
bound is given by the difference between the volume fraction at the surface and 
the center of the film. Figure 3 demonstrates that the theoretical results for 
A<&a lie nicely within the bounds provided by the simulation data. However, 
quantitative agreement is only reached if the range of monomer interactions 
has been accounted for correctly, i.e. k 2 may not be neglected. 

The total amount of A segregated at the surface is given by the total excess 
Z\ = J dzp(z)(&A(z) — &A,b)- As shown in Figure 4, it can be described by a 
simple Boltzmann dependence 

Z* A oc exp(-xiVa) (25) 

over a wide range of xN- In this regime, one can describe the segregation 
process as adsorption of single, relatively undistorted A chains whose energy 
is reduced by xN(l — °0 & t the surface, compared to the bulk. The quantity 
(1 — a) thus gives the quota of neighbor monomers, which a chain is missing 
at the surface. In the present model, one obtains a = 0.77. Close to the 
critical point, interactions between different A chains at the surface become 
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important and our simple picture breaks down - the surface excess diverges. 
At very high values of xN, beyond x^^14, the chains are more and more 
squeezed at the surface, the number of missing neighbors increases, and the 
segregation is slightly higher than eqn (|25| ) would predict. 

The surface segregation in the one phase region is shown in Figure 5, as a 
function of bulk volume fraction $^,6 at fixed Flory Huggins parameter x- For 
comparison, segregation profiles for mixtures of polymers with different chain 
lengths (N B = 2N A = 48 or N A = 2N B = 48 as opposed to N B = N A = 32) 
were also calculated. In such asymmetric mixtures, shorter chains segregate to 
the surface when there are no interactions, at x = 0, for entropic reasons (thin 
line). At x'ZXc however, the segregation is mainly driven by energetic effects: 
As the bulk volume fraction $^ approaches zero, polymers A segregate to the 
surface to the same extent regardless of whether they have longer or shorter 
chain lengths than polymers B. The effects of segregation due to a difference 
in chain lengths, and segregation due to missing neighbor effects, are therefore 
not additive. This becomes most obvious when looking at the surface excess of 
component A (Figure 5b). On increasing the fraction of chain lengths N A /N B 
while keeping Xc = (1/v^Va + 1/ V 'N B ) 2 /2 approximately constant, the critical 
point moves towards the A-poor region of the phase diagram (cf eqn (f24|)). 
As a consequence, the excess of A-polymers at the surface is even higher for 
N A > N B than for N B > N A . Hence the entropic mechanism which promotes 
segregation of the shorter chains is largely suppressed - the dominant effect of 
chain length asymmetry is the shift of the critical point. 



3.2 Structure of the Profiles 

We turn to analyze the segregation profiles in more detail. Far from the critical 
point, they can be expected to reflect single chain behavior. Figure 4 suggests 
that one may be able to describe a segregation profile in terms of a weakly 
adsorbed layer of the minority component. In order to test this picture, we 
calculate for comparison the density profile of a single Gaussian polymer in 
semi-infinite space, which touches the surface at least once. This can be done 
analytically, as shown in the appendix. Far from the surface, the density pro- 
file decays asymptotically like p(z) oc 1/z 3 exp(—z 2 /(AR 2 )). The logarithmic 
derivative is thus given by 

R g dl ° gp{z) = -™3L - (26) 
9 dz z 2R g K ' 

Figure 6 shows the logarithmic derivative of segregation profiles p A (z) in the 
two-phase region at different values of x, compared to the prediction of eqn 
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(|26| ) (dashed line). Far from the critical point, our simple picture describes 
the asymptotic behavior surprisingly well. Close to the critical point, the log- 
arithmic derivative of Pa{ z ) approaches a constant at large z: The segregation 
profile decays exponentially with a decay length £. This is characteristic for 
a system close to a critical point. The interactions between chains become 
important, and collective effects determine the structure of the profiles. In 
polymeric systems, the characteristic length scale £ for collective concentra- 
tion fluctuations competes with the gyration radius R g , which describes the 
correlations produced by chain connectivity - the correlations of monomers 
belonging to the same chain. The result is a crossover between two types of 
qualitatively different asymptotic behavior: Far from the critical point, intra- 
chain correlations dominate, close by, collective concentration correlations do. 
The correlation length £ diverges at the critical point. 

We can again attempt to compare these results with the simulations of 



Ref. |]18[ . Simulations where performed at values xN = 3 — 4, where the 
asymptotic behavior of the profiles is still governed by the length scale £. At 
film thicknesses of less than 3R g , however, the asymptotic behavior is never 
reached. This may explain why the decay length obtained by Rouault et al, 
£ ~ 0.35R g at xN — 3.4, is only roughly half the value predicted by the self 
consistent field theory, £ « 0.6R g . 

The decay length £ as a function of 1/xiV is shown in Figure 7, and com- 
pared to half the width of the bulk interface between the two coexisting phases, 
which has been calculated in Ref. P5| . In a simple system, both lengths are 
identical, i.e. by the bulk correlation length. Here, this is only true very close 
to the critical point; further away, £ gets too large relative to the interfacial 
width. This is due to the fact that the wings of the segregation profiles are 
governed by the chain end distribution, as we shall see in the next section. 
The interfacial width, on the other hand, describes the separation of all A and 
B monomers. Chain end distributions tend to be broader than total monomer 
distributions for entropic reasons. In the wings of the interfacial profile, we 
expect to find a decay length which is again determined by the chain end 
distributions and identical to £. This is indeed the case (not shown). 

Figure 7 also shows data for another length which characterizes the segre- 
gation profiles, the profile width. It is conveniently defined through the excess 
distribution Z A (z) = p(z)(Q>A.(z) — &A,b)/Pb, e.g. as its first moment or as its 
standard deviation. 



Jdz zZ A (z) 



W {z) = W {SZ 2) 



Jdz z 2 Z A (z) _ Jdz zZ A (z) ^ 

\ ^A 

At the critical point, the width diverges with the correlation length £. Further 
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away, it reaches a plateau at W ~ R g f a/3, which is about the extention of an 
undistorted chain in one dimension. At large values of x-W, however, the width 
decreases further and approaches zero. The strong repulsion between chains A 
and B cause the segregated chains to get increasingly distorted and squeezed 
towards the surface. The turning point of the curve is at l/%iV ~ 0.07 or 
\N « 14, consistent with what was already deduced from the x-^~dependence 
of the surface excess Z* A in Figure 4. 



3.3 Chain end distributions and Chain orientations 

Figure 8 compares the distribution of chain ends with density profiles close 
to the surface at bulk two phase coexistence. For both the minority and the 
majority component, the fraction of chain ends relative to the total density is 
augmented right at the surface, and reduced in a region underneath. Such an 
effect is already known from pure systems [O, [TH]. In the minority component, 



the depletion zone is followed by another enrichment zone in the wings of the 
decaying segregation profile - polymers of the minority component rather stick 
out chain ends into the majority phase than whole loops. Thus the segregated 
layer of polymers A is skirted by regions of enhancement of chain ends A. 
The total chain end distribution profile, however, only has one enrichment 
zone right at the surface. One can probably expect a second zone at higher 
segregation, beyond the wetting transition, since chain ends are known to 



enrich at the A/B interface |43|, In our simple picture of an adsorbed 
polymer layer, the chain end distribution decays asymptotically like p(z) oc 
1/z exp(— z 2 /(4Rg)) (see appendix), i.e. the logarithmic derivative follows 

d\ogp end (z) = _R g _ j_ 
9 dz z 2R g 1 ' 

This is compared to the logarithmic derivative of chain end distributions of 
the segregating component at different values of \ in Figure 9. The agreement 
for high values of x~N is not as good as for the total density profiles, but still 
remarkable. Close to the critical point, the profiles show asymptotically again 
exponential decay with the decay length £. 

The orientation of whole chains can be studied by solving the diffusion 
equation (j|) in the previously determined self consistent field Wi with initial 
condition Qi(f, 0) = S(f—r ). The statistical weight of a polymer with end-to- 
end vector R e = r*i— r and center at R c = (fi+ro)/2 is then given by Qi(fi, 0). 
Since the self consistent field W{ only varies in the direction perpendicular 
to the interface, z, the components of the end-to-end vector parallel to the 
interface have the same distribution at the surface than in the bulk. The 
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presence of a surface only affects the ^-component R Zt6 . Chains get therefore 
oriented as a result of stretching or squashing in the z direction, the theory 
does not allow for the possibility of chain orientation without compression, as 
has been observed in computer simulations Wffl. Profiles of {R e , z 2 )(z), where 
the averages have been performed for polymers with center at z, are shown in 
Figure 10. The bulk value of this quantity is simply 2R g 2 with the gyration 
radius R g . Near the surface, it is greatly reduced - it drops down to zero at 
z = and remains smaller than the bulk value over a distance of more than 
one gyration radius from the surface. Hence polymers tend to align parallel to 
the surface. This region is followed by a zone of perpendicular alignment at 
the distance 2R g from the surface (Figure 10). Both the parallel alignment [|15H 
and the subsequent "overshoot" of perpendicular alignment [|16| have been seen 
in Monte Carlo simulations of on one component polymer melts. The latter 
effect is weak for chains of the majority component, but can be quite marked 
for the minority component, deep in the two-phase region. The extent of 
parallel alignment, on the other hand, does not depend on whether a chain 
belongs to the minority or the majority component, or on the strength of the 
interaction \. Thus we find that the total structure of the blend does not 
change with increasing x, although the segregated minority component chains 
get more and more distorted, as indicated by Fig. 7. They just accumulate 
closer to the surface, where the distortion is higher. 



4 Summary 

A self consistent field theory for compressible binary polymer mixtures has 
been developed, which allows for a consistent treatment of surfaces, including 
effects of missing neighbors. It has been shown to be remarkably successful 
in reproducing Monte Carlo results for the segregation at a neutral surface in 
the bulk coexistence region, without resorting to any fitting parameter. On 
this base, several aspects of the surface structure could be discussed, like the 
role of chain end distribution and chain orientations, the effect of chain length 
asymmetry, and the interplay of the different length scales which govern the 
system: the correlation length of collective concentration fluctuations and the 
single chain correlation length, i.e. the gyration radius. Hence we have seen 
that the formalism developed in this work is suited for the investigation of 
surfaces of polymer mixtures. So far it has only been applied to systems with 
perfectly symmetric interactions (where V^{r) oc {Uaa(j) — Ubb(j)) (eqn (|16|)) 
vanishes), a rather unlikely situation in real systems. Usually the asymmetry of 
the interactions introduces additional surface potentials h(z) (cf. eqn (pT), i.e. 
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preferential attraction of one component to the surface, After taking those into 
account, the theory should be similarly successful in describing experimental 
results. 

I wish to thank K. Binder for suggesting the problem and for useful com- 
ments on the manuscript, and L. Klushin and O. Borisov for helpful discussions. 



Appendix 



We consider a Gaussian Polymer of length N in the half space with end 
monomers located at distances zq and Z\ from the surface, which touches the 
surface with the K\h. monomer (Figure 11). The Greens function for propa- 
gation in the half space can be obtained with the method of mirror images 

(z-z') 2 (z + z') 2 

G(M, z; M', z') oc expf-^ 7777.7) - exp(- 2/ \, A, to ), (29) 

where b is the Kuhn length, and z, z' are the locations of monomer M and M' 



[48]. Hence the distribution profile of monomer L is given by 

D L (z\0,z ;K,e^0 + ;N,z 1 ) oc (30) 



lim 



G(0,z ;L,z)G(L,z;K,e)G(K,e;N,z 1 ) : K > L 
0+ \ G{0, z ; K, e)G(K, e; L, z)G{L, z; N, z\) : L > K ' 

In order to get the total density of the Lth monomer, one has to integrate over 
the free parameters Zq, z\ and K. The result is 

p L (z) oc gi(az)VlVl - I (31) 

with 



gi(x) = -x 



$(^)Ei(-^-)— = + $(^^) E i(-^)^ 



where / = L/N and a = l/(2R g ) with the gyration radius R 2 = Nb 2 /6. The 
error function $ is defined as $(x) = ^= J a: exp(— t 2 ) dt, and the exponential 
integral Ei as — Ei(— x) = exp(—t)/tdt. Note that the normalization factor 
Af(Vl) = / °° dxgi(x) is finite and nonzero for all values of / between zero and 
one. The total density profile is given by 

p{z) = f 1 dl gi (az)/M(Vl). (32) 



o 
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The asymptotic behavior of the profiles can be obtained from the asymptotic 
behavior of gi(x): 



9i(x) ™ \[ exp(- I ^ 7 ) v / r^T + exp(-^)v / / ]. (33) 
The result for the chain end distribution follows immediately. 

Po(z) = pi{z) oc - exp(-(az) 2 ) (34) 

z 

The asymptotic behavior of p(z) is slightly more difficult to derive. Wihout 
calculating M~ 1 {y) explicitly, we expand it in powers of y: M~ x {y) = J2k °ky k - 
Using asymptotic series representations of the error function and the exponen- 
tial integral, one can show 

I * ^ ~ + (35) 

with the Gamma function Y and TZn < x~ 2( - N+1 ^T(N). With that one obtains 

o oo r i 

az k=Q Jo 

M 3 fc ^o ( } r((fc + i)/2) ^ 



i 



2c b e-<"> — (36) 
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Figure Captions 



Figure 1: Mean field phase diagram of the symmetric polymer blend with 
Na = Nb = N in the plane of variables 1 / (xN) and the volume fraction 
$a of A monomers. Inside the binodal curve (solid line), the mixture 
phase separates into two macroscopic phases (1) (left) and (2) (right). 
Also shown is the surface volume fraction of A monomers at a free sur- 
face in each of the two coexisting phases. Near the critical point, 

the concentration differences ($^ — and {&a\ — &a\) vanish with 
different critical exponents: — ^a) ^ (Xc — xY with (3 = 1/2 and 

(&a\ ~ ^a\) (Xc — xY 1 with f3\ = 1 in mean field theory. Thus the 
broken curves have a cusp at the critical point. 

Figure 2: Volume fraction profile of the minority component monomers at 
x/Xc — 1-696 in the symmetric polymer mixture. Dashed line shows 
volume fraction profile for ft2 = 0, where monomers are assumed to have 
pure contact interactions with each other. Inset shows the total density 
profile, which was identical for all considered values of x- 

Figure 3: Difference between surface and bulk volume fraction of the 

minority component monomers in the symmetric polymer mixture, at two 
phase coexistence, vs. l/(%./V). Data points show results of the simula- 
tions of Rouault et al (flTB|) : The difference between surface and bulk 
volume fraction (circles), and the difference between the surface volume 
fraction and the volume fraction at the center of simulated slab (squares). 
The arrow indicates the value of 1/ (xN) corresponding to the profiles in 
Figs 2 and 8. 

Figure 4: Surface excess Z* A of the minority component at two phase coexis- 
tence in the symmetric polymer mixture, in units of PbR g , vs \N- dashed 
line shows the function 0.3 exp(— 0.77xiV). 

Figure 5: (a) Difference between surface and bulk volume fraction of A 
monomers and (b) surface excess Z* A in units of ptR g , in the one phase 
region, as a function of the bulk volume fraction Considered are symmet- 
ric mixtures with Na = Nb (thick solid line) and asymmetric mixtures, 
where A polymer chains are either longer (long dashed line) or shorter 
(short dashed line) than B polymer chains. Thin line indicates results for 
asymmetric mixtures with the interactions turned out or x — 0- 

Figure 6: Logarithmic derivative of density profiles for the minority compo- 
nent in the symmetrical mixture, in units of R~ l , at two phase coexistence. 
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Values of x/Xc are, from top to bottom: 

X / Xc = {1.102, 1.271, 1.695, 2.542, 3.390, 4.237, 5.085, 5.932, 6.780, 7.627, 8.475}. 
Dashed line shows the function —z/(2R g ) — (3R g )/z for comparison (see 
text for explanation). 

Figure 7: Characteristic lengths of the segregation profiles in units of the 
gyration radius, vs l/(%iV), in the symmetrical mixture at two phase 
coexistence. Solid line and long dashed line show the width of the profiles, 
defined in two different ways. Short dashed line gives asymptotic decay 
length £, dotted line gives half the interfacial width between the two 
coexisting phases in the bulk. Dashed dotted line indicates the extension 
of one undistorted chain in the bulk. Triangles mark the values of 1/ (%iV) , 
corresponding to the profiles shown in Figs 6 and 9. 

Figure 8: Distribution of chain end monomers of the minority component 
PA,end i n units of (2/0{,/jV) in the symmetric polymer mixture at two 
phase coexistence, x/Xc — 1-696. Dashed line indicates distribution of 
all monomers for comparison. Inset shows the total chain end distribu- 
tion (solid line), compared to the total density profile (dashed line). 

Figure 9: Logarithmic derivative of chain end distributions for the minority 
component in the symmetrical mixture, in units of Rj 1 , at two phase 
coexistence. Values of x/Xc are as in Figure 6, from top to bottom: 
X/Xc = {1.102, 1.271, 1.695, 2.542, 3.390, 4.237, 5.085, 5.932, 6.780, 7.627, 8.475}. 
Dashed line shows the function —z/{2R g ) — R g /z 

Figure 10: Profiles of the average squared ^-component of the end-to-end 
vector {R 2 ez ) for chains centered at z/R g , vs z/R g . Considered are chains 
of the minority component (dashed lines) and the majority component 
(solid line) in a symmetrical mixture at two phase coexistence, different 
values of x- The results for chains of the majority component do not 
depend on x- 

Figure 11: Cartoon of a single weakly adsorbed chain in semi-infinite space 
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